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Convective  heat  transfer  is  important  for  maintaining  life  and  reliability  of  the  main 
combustion  chamber  in  liquid  rocket  engines.  Enhancing  the  regenerative  cooling  circuit’s 
ability  to  maintain  wall  temperature  with  minimal  pressure  losses  improves  engine  life  and 
has  the  potential  to  reduce  the  turbopump  output  pressure.  High  aspect  ratio  coolant 
channels  have  been  shown  to  reduce  wall  temperature  with  little  pressure  increase. 
Additionally,  rocket  engine  cooling  circuits  are  curved  and  the  mechanisms  responsible  for 
curvature  enhancement  remain  unclear  for  tall,  narrow  passages  and  appropriate  predictive 
tools  have  not  been  developed.  The  objective  of  this  research  is  to  study  the  mechanisms 
responsible  for  improved  cooling  performance  and  to  develop  heat  transfer  correlations  for 
use  in  high  aspect  ratio  coolant  channels  having  conducting  side  walls  and  subjected  to 
asymmetric  heating  on  the  concave  side,  thus  modeling  a  channel  in  a  slotted  liner.  Tests  are 
conducted  with  water  in  a  single,  curved  channel  subjected  to  heat  flux  on  the  outer  radius 
wall  and  in  a  straight  channel  with  heat  input  to  an  analogous  wall.  Analyses  are  done  in 
three  areas.  First,  correlations  are  developed  as  functions  of  Reynolds  number,  Prandtl 
number,  ratio  of  viscosities,  aspect  ratio,  axial  location,  and  curvature  based  on 
experimental  data.  Analysis  shows  that  approximately  96%  of  curved  channel  data  are 
predicted  to  within  +10%  by  the  resulting  Nusselt  number  correlations.  The  average 
percent  difference  between  the  data  and  predicted  value  is  approximately  3%  of  the 
predicted  value.  This  indicates  that  the  correlations  properly  capture  the  trends  and  yield 
useful  predictions  of  heat  transfer.  Secondly,  published  correlations  are  evaluated  at  each 
angular  location  capture  current  data  trends  well.  Lastly,  analysis  examines  effect  of 
curvature  on  pressure  gradients.  Computational  fluid  dynamics  is  incorporated  to  provide 
insight  into  mechanisms  responsible  for  pressure  losses.  Radial  pressure  gradients  start 
forming  approximately  11  hydraulic  diameters  upstream  of  the  start  of  curvature  and 
continue  approximately  18  hydraulic  diameters  downstream  of  the  end  of  curvature.  The 
data  found  in  this  investigation  is  of  good  quality  and  will  be  useful  when  moving  to  subscale 
testing. 
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Nomenclature 

AR  =  Aspect  Ratio 

C,  m,  n,  q,  r,  s,  u-  Unknown  constants  in  the  Nusselt  Number  correlations 

dc  =  Diameter  of  Curvature 

De  =  Dean  Number 

Dh  =  Hydraulic  Diameter 

dP/dr  =  Radial  Pressure  Gradient 

ge  =  Earth’s  Gravitational  Constant 

g  *  =  Centripetal  Acceleration  (non-dimensionalized) 

H  =  Channel  Height 

h  =  Average  Heat  Transfer  Coefficient 

k  =  Thermal  Conductivity 

L  =  Length 

Nuc  =  Nusselt  Number  for  curved  channel  data 
Nus  =  Nusselt  Number  for  straight  channel  data 
q  ”  =  Heat  Flux 

R0  =  Radius  to  the  bottom  of  the  channel 

Re  =  Reynolds  Number 

Pr  =  Prandtl  Number 

5  =  Arc  Length 

Tb  =  Bulk  Temperature 

=  Inlet  Temperature 

U  =  Bulk  Averaged  Velocity 

v  =  Axial  Distance 

W  =  Channel  Width 

6  =  Angular  Distance 

(pcur  =  Curvature  term  in  Nusselt  Number  correlations 


I.  Introduction 

Proper  design  of  the  cooling  circuit  of  a  liquid  rocket  engine  is  critical  for  maintaining  life  and  reliability  of  the 
combustion  chamber.  The  most  challenging  location  for  accomplishing  this  design  is  the  throat  region  due  to 
the  elevated  hot-gas-side  heat  flux,  geometrical  constraints,  and  curvature  of  the  coolant  channel  flow  path. 
Fortunately,  the  curvature  provides  some  enhancement  to  the  coolant-side  convection  coefficient  over  that  which 
would  be  present  in  a  straight  channel  (Seban  and  McLaughlin,  1963;  McCormack,  Welker,  Kelleher,  1969;  Sturgis 
and  Mudawar,  1999).  However,  it  is  unclear  whether  the  mechanisms  responsible  for  this  enhancement  exist  when 
the  coolant  channel  assumes  a  high  aspect  ratio  (i.e.  tall  and  narrow  cross-section).  Wadel  (1998)  showed  that  the 
hot-gas  wall  temperature  could  be  lowered  with  only  small  increases  in  pressure  losses  when  examining  high  aspect 
ratio  coolant  channels  (HARCC)  using  an  analytic  method  which  was  anchored  by  hot-fire  data.  All  HARCC 
designs  considered  lowered  hot-gas  wall  temperature  by  at  least  8%  with  as  little  as  2%  increase  in  pressure  loss. 
Wadel  concluded  that  HARCC  has  an  overall  benefit,  regardless  of  coolant  channel  configuration  studied. 

Implementation  of  HARCC  is  becoming  more  feasible  since  machining  techniques  such  as  milling  are  now  able 
to  fabricate  deep,  narrow  grooves  of  varying  cross-section  on  contoured  surfaces.  Therefore,  the  issue  from  a  design 
standpoint  is  to  take  advantage  of  both  the  curvature  enhancement  and  high  aspect  ratio  geometry  of  the  coolant 
channels  with  adequate  predictive  tools  and  with  proper  insight  into  the  physics  of  the  interaction. 

Relatively  few  studies  have  focused  on  quantifying  the  effects  of  high  aspect  ratio  coolant  channels  and  even 
fewer  on  the  effects  of  HARCC  on  curvature  enhancement.  Also,  the  effects  of  HARCC  on  curvature  enhancement 
mechanisms  are  not  well  understood.  Therefore,  the  combination  of  lack  of  data  and  unclear  understanding  of 
HARCC  presents  a  design  risk  when  using  HARCC.  Available  heat  transfer  correlations  were  developed  for  tubular 
passages  however  they  fail  to  capture  the  particular  characteristics  of  high  aspect  ratio  channels  in  a  rocket  engine 
cooling  circuit;  namely,  rectangular  passage,  one-sided  heating,  temperature  gradient  in  side  wall  and  flow 
curvature.  It  is  important  that  the  thermal  designer  understand  how  the  enhancement  is  influenced  by  the  channel 
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geometry  and  be  able  to  predict  this  for  different  conditions.  Since  the  engine  throat  represents  the  most  severe 
thermal  environment  in  the  chamber,  the  consequences  of  inadequate  design  tools  could  caused  significant 
reductions  in  life  and  reliability,  if  not  catastrophic  failure. 

The  current  research  uses  a  synergistic  experimental  and  computational  approach  to  study  curved  and  straight, 
asymmetrically-heated,  high  aspect  ratio  coolant  channels.  The  overall  effort  focuses  on  three  aspects: 
understanding  the  internal  mechanisms;  generating  new  heat  transfer  correlations;  and  validating  CFD  results  in 
order  to  improve  the  predictive  tools.  Moreover,  curved  channel  data  are  compared  to  the  baseline  straight  channel 
experimental  data  permitting  the  quantification  of  the  curvature  enhancement  factor  in  HARCC.  Curved  channel 
experimental  data  are  incorporated  in  this  paper  from  only  aspect  ratios  4  and  10,  though  other  tests  are  ongoing. 
This  paper  focuses  on  generating  heat  transfer  correlations  for  curved  channel  data,  comparing  present  correlations 
to  straight  channel  correlations,  assessing  several  published  curvature  terms,  and  investigating  curvature  effects  on 
pressure  losses. 


II.  Experimental  Methods 

The  High  Aspect  Ratio  Coolant  Channel  (HARCC)  Water  Flow  Loop  at  the  Air  Force  Research  Laboratory  is 
used  for  testing  single  straight  and  curved  channels  subjected  to  asymmetric  heating.  Open  to  atmosphere,  the 
recirculation  loop  permits  continuous  testing  while  maintaining  constant  inlet  conditions  (temperature,  pressure,  and 
flow  rate)  to  the  test  section. 

A  typical  cross-section  of  the  curved  copper  test  channel  is  shown  on  in  Figure  1  for  aspect  ratio  AR  =  4. 
Straight  test  channel  has  identical  cross-section.  The  channels  are  formed  by  cutting  a  slot  of  width,  W  =  6.35  mm 
into  a  single  piece  of  102  copper  alloy.  A  T-insert  is  placed  in  the  slot  to  form  a  rectangular  flow  channel.  The 
length  (L  =  381mm)  along  the  bottom  of  the  rectangular  flow  channel  is  the  same  for  both  the  curved  and  straight 
channels.  For  example,  the  arc  length  along  the  bottom  curved  channel  wall,  S  =  Rqk  =  L  equals  the  heated  length 
of  the  straight  channel,  where  R0  =  121mm.  For  the  current  test  matrix,  five  such  T-inserts  are  used  to  create  aspect 
ratios  of  AR  =  1,  2,  4,  8,  and  10  which  correspond  to  heights  of  H-  6.35,  12.7,  25.4,  50.8,  and  63.5  mm.  Each  of  the 
side  walls  is  half  the  channel  width  and  the  copper  is  completely  enclosed  in  insulation.  The  result  is  a  rectangular 
flow  channel  with  conducting  side  walls  that  models  a  single  channel  in  a  regenerative  cooling  circuit  whose  land 
width  is  equal  to  the  channel  width. 

The  test  section  consists  of  two  portions:  an  aluminum  entry  channel  and  then  a  copper  test  channel  which  can  be 
straight  or  curved.  For  more  details  about  the  straight  channel  see  the  technical  memos  by  Nathman  and  Sturgis  in 
2008.  Figure  2  shows  the  entry  channel  with  the  curved  copper  channel  installed  in  the  flow  loop  ready  for  testing 
with  the  exception  of  insulation  which  has  been  removed  for  clarity  in  illustrating  the  test  rig.  The  unheated  entry 
channel  provides  a  fully  developed  turbulent  flow  into  the  heated  copper  test  section  with  a  length  of  508  mm 
corresponding  to  80  hydraulic  diameters  ( Dh )  in  length  for  AR  =  1  and  44 Dh  for  AR  =10.  In  any  given  test 
configuration,  both  channels  have  the  same  internal  passage  geometry. 

Heat  is  applied  at  the  outer  radius  by  41  cartridge  heaters  embedded  in  the  base  of  the  copper  block  (Figs.  1  &  2) 
connected  to  a  variable  voltage  source.  Due  to  facility  power  limitations,  heat  fluxes  up  to  only  130  W/cm2  (0.8 
Btu/in2s)  could  be  reached  in  the  curved  channel  whereas  180  W/cm2  (1.1  Btu/in2s)  is  possible  in  the  straight 
channel.  Future  upgrades  will  mitigate  this  problem.  Heat  flow  to  the  channel  is  modeled  as  one-dimensional  in  the 
constant  area  cross-section  since  the  copper  is  well  insulated. 

Thermocouples  embedded  in  the  cross-section,  shown  in  Fig.  1,  permit  the  calculation  of  heat  flux  and 
convection  coefficient.  Three  thermocouples  in  the  thick  neck  section  determine  the  one -dimensional  heat  flux 
acting  on  the  channels.  Fin  thermocouples  indicate  heat  flow  up  the  side  walls  and  into  the  fluid  and  are  used  in  the 
estimation  of  an  average  convection  coefficient  for  the  cross-section.  The  array  of  thermocouples  exists  at  five 
angular  locations  in  the  channel  which  are  indicated  in  Fig.  1  and  visible  in  Fig.  2.  This  permits  calculation  of  heat 
flux  and  convection  coefficient  along  the  flow.  Subsequent  discussions  will  focus  on  three  locations  -  channel  inlet 

( 6  =  9°),  mid-channel  ( 0  =  90°),  and  channel  exit  {0  =  171°)  -  which  are  deemed  sufficient  to  capture  the  trends. 

These  angular  locations  correspond  to  axial  locations  (x  =  19.1,  191,  and  364  mm)  on  the  straight  channel.  Pressure 
taps  are  located  near  the  inlet  and  exit:  6  =  15°  and  165°  for  the  curved  channel  or  L  =  31.75  mm  and  349.25  mm  for 
the  straight  channel.  The  pressure  tap  locations  allow  for  measurement  of  the  radial  pressure  gradient,  dP/dr,  for  AR 
=  8  and  10.  For  lower  aspect  ratios  the  T-insert  interferes  with  the  pressure  taps  at  the  top  of  the  channel,  so  the 
pressure  gradient  cannot  be  measured. 
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Figure  1.  Curved  Copper  test  channel  showing  thermocouples  at  five  angular  locations  and  the  four 

pressure  taps 
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Aluminum  Curved  Copper 

Entry  Length  Test  Channel 

Figure  2.  Photo  of  aluminum  entry  channel  and  curved  copper  test  channel  connected  to  water  flow 

loop 
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A  typical  test  series  begins  by  configuring  the  test  channel  (and  entry  channel)  with  a  certain  aspect  ratio  and 
inserting  it  into  the  flow  loop.  After  establishing  inlet  pressure,  temperature,  and  flow  rate,  the  heat  input  is 
initiated.  Data  is  recorded  at  successively  higher  levels  of  heat  flux,  in  each  case  waiting  for  steady  state  conditions 
to  prevail.  The  test  concludes  at  the  highest  heat  fluxes  or  when  the  channel  wall  temperature  approaches  100  °C. 

A  new  test  is  initiated  by  adjusting  the  flow  rate  to  achieve  a  different  average  velocity.  Nominally,  velocities  of 
U  =  5,  10,  and  15  m/s  are  achieved.  These  tests  are  accomplished  for  two  inlet  fluid  temperatures,  Tin  =  20  and  35 
°C.  This  series  of  tests  is  repeated  for  each  of  the  five  aspect  ratios  examined,  with  several  extra  tests  to  assess 
repeatability.  Therefore,  the  experimental  heat  transfer  data  is  collected  for  ranges  of  axial  location  (1.7  <  x/Dh  < 
57),  heat  flux  (10  <  q”  <  130  W/cm2),  Reynolds  number  (50,000  <  ReDh<  300,000),  fluid  temperature  (Tin  =  20, 
35  °C)  and  aspect  ratio  (AR  =  4  and  10).  The  inlet  pressure  is  held  constant  at  2.1  bars  for  all  tests.  Measurement 
uncertainties  entail:  0.7  °C  for  Tb;  0.0086  bar  (0.125  psi)  for  pressure;  3.5%  max  for  U\  30%  for  low  q”  (15  W/cm2) 

decreasing  to  13%  or  less  for  q”  >  40  W/cm2  which  represents  most  data;  h  and  Nu  will  be  slightly  more  than  q” 
due  to  calculation  assumptions;  Re  uncertainty  is  similar  to  U\  dP/dr  uncertainty  is  similar  to  pressure. 

Several  potential  parameters  governing  curved  flow  have  been  identified.  The  first  is  purely  a  geometric  term  - 
the  ratio  of  diameter  of  curvature  to  diameter  of  flow  passage.  It  is  less  clearly  defined  for  a  rectangular  passage 
since  the  diameter  of  curvature  ( dc )  could  be  interpreted  as  twice  the  distance  from  the  center  of  curvature  to  the 
center  of  the  cross-section  or  to  the  heated  wall.  The  difference  can  be  significant  with  high  aspect  ratio  coolant 
channels.  For  this  study  dc  =  2R0  is  used.  The  diameter  of  the  flow  passage  is  usually  taken  as  the  hydraulic 
diameter.  The  second  is  the  Dean  number  which  represents  the  ratio  of  the  square  root  of  the  product  of  inertial  and 
centrifugal  (pressure)  forces  to  the  viscous  forces  (Berger  et  al.,  1983).  A  third  term  is  the  centripetal  acceleration 
experienced  by  the  fluid  non-dimensionalized  by  earth’s  gravitational  constant,  ge.  This  also  is  subject  to  selection 
of  the  curvature  dimension,  but  for  the  present  work  will  be  defined  as 


g*  = 


(i) 


where  R0  is  the  radius  at  the  channel  outer  wall.  The  ranges  of  these  parameters  are:  21  <  2 RJDh  <  24;  6,000  <  De 
<  66,000;  and  21  <  g*<  190. 

It  should  be  noted  that  the  present  test  conditions  do  not  represent  normal  operating  environments  for  rocket 
engine  coolant  channels.  For  example,  representative  Reynolds  number  and  heat  flux  values  for  rocket  engine 
cooling  circuits  exceeds  lxlO6  and  8000  W/cm2  (~  50  Btu/in2-s),  respectively.  However,  the  objective  is  to 
understand  trends  and  perform  tests  at  a  more  feasible  level.  Computational  efforts  will  be  used  to  extend  the 
insights  and  correlations  developed  through  these  tests  to  more  realistic  rocket  conditions.  Additionally,  another 
research  program  is  under  way  to  examine  HARCC  by  means  of  a  small  combustion  chamber  which  can  achieve 
higher  Reynolds  numbers  and  heat  fluxes.  These  tests  will  complement  the  present  investigation  and  validate  the 
computational  results.  The  test  section  is  instrumented  such  that  validation  and  appropriate  analysis  can  be  done. 

Data  reduction  proceeds  from  measured  temperatures,  flow  rates,  pressures,  and  key  assumptions.  Average 
velocities  are  calculated  from  measured  flow  rate  and  channel  cross-sectional  area.  Fluid  inlet  temperature  is 
measured  by  a  thermocouple  inserted  into  the  flow.  Bulk  temperature  (Tb)  at  each  of  the  axial  instrumentation 
locations  is  calculated  based  on  known  flow  rate,  fluxes  up  to  the  point  of  interest,  and  the  assumption  of  a  well- 
mixed  flow.  Water  properties  such  as  density,  viscosity,  thermal  conductivity,  and  Prandtl  number  are  considered 
functions  of  only  bulk  temperature.  Applied  heat  flux  ( q  ”)  derived  from  the  copper  temperature  gradient  in  the  neck 
region  (indicated  by  thermocouples)  and  assumption  of  constant  copper  conductivity  (390  W/m-K). 

The  convection  coefficient  is  derived  from  an  energy  balance,  measured  temperatures  and  certain  assumptions. 
The  applied  heat  load  (in  neck  region)  is  equated  to  the  heat  transferred  to  the  fluid  by  convection  plus  that 
conducted  up  the  fin  beyond  the  flow  passage.  Estimations  indicate  that  losses  through  the  side  insulation  are 
negligible.  Convection  to  the  fluid  is  estimated  by  integration  along  the  channel  surface  using  measured  fin  (wall) 
temperatures  and  assumption  of  a  bulk  fluid  temperature.  Conduction  up  the  fin  beyond  the  flow  passage  is 
estimated  by  the  temperature  gradient  in  the  fin  at  the  top  of  the  flow  passage  and  is  negligible.  The  energy  balance 

yields  an  average  convection  coefficient,  h  ,  for  the  cross-section  where  it  is  understood  that  h  pertains  to  only  three 
walls.  A  convection  coefficient  is  determined  for  each  flux  value  recorded  at  three  axial  locations  -  near  the  inlet, 
midway  along  the  channel,  and  near  the  exit.  Corresponding  Nusselt  numbers  are  calculated  by  multiplying  the 
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convection  coefficient  by  hydraulic  diameter  and  dividing  by  the  water  conductivity.  For  detailed  descriptions  of 
the  experimental  setup  and  data  reduction  refer  to  technical  memos  by  Sturgis  and  Nathman  in  2008. 

III.  Results  and  Discussion 

This  initial  examination  of  AR  =  4  and  10  experimental  heat  transfer  data  is  used  to  develop  Nusselt  number 
correlations  for  curved  flow,  compare  curved  channel  results  to  straight  channel  results,  and  evaluate  published 
curvature  terms.  Secondly,  pressure  gradient  trends  in  curved  flow  are  analyzed. 

A.  Correlation  Analysis 

The  objective  of  developing  Nusselt  number  correlations  is  to  create  a  design  and  analysis  tool  for  the  thermal 
engineer.  A  correlation’s  usefulness  depends  on  its  ability  to  accurately  predict  the  heat  transfer  coefficient  for  the 
conditions  being  examined.  These  conditions  include  channel  geometry,  heating  arrangement,  fluid  characteristics 
and  axial  location. 

In  general,  the  Nusselt  number  is  correlated  with  Reynolds  and  Prandtl  numbers  following  the  basic  form 
proposed  by  Dittus  and  Boelter  (1930).  In  this  expression, 

Nu  =  C  RemPrn  (2) 

the  unknowns  C,  m,  and  n  are  determined  from  a  least  squares  fit  of  the  experimental  data.  The  current  work  takes 
the  initial  form  and  builds  on  it  to  eventually  create  the  expanded  expression, 

( Y  (  d  y 

Nu  =  C  Rem  Prn  ^  ARq  1  +  -^—  cpucur  (3) 

/  7-^  /  CUT  v  ' 

VMW)  {  x/D  ) 

where  unknowns  C,  m,  n,  p,  q,  r,  and  s  are  determined  from  a  least- squares  best-fit  to  the  experimental  data. 
Equation  3  includes  Reynolds  number  (Re),  Prandtl  number  ( Pr ),  ratio  of  water  viscosity  evaluated  at  bulk 
temperature  and  channel  base  wall  temperature  (/v^w),  aspect  ratio  (AR),  thermal  flow  development  (1  +B/( x/D)), 
and  the  curvature  parameter  (q>cur).  Less  critical  is  the  value  of  B  which  is  chosen  in  order  to  avoid  large  values  of 
the  bracketed  term  at  the  channel  inlet  based  on  the  present  x/Dh  data.  In  this  paper  only  results  from  AR  =  4  and  10 
will  be  presented.  Since  correlations  will  be  developed  for  each  aspect  ratio,  the  AR  term  is  omitted  from  this  best- 
fit  analysis  due  to  its  singular  value.  Due  to  delays  in  analysis,  the  value  of  exponent  u  and  the  curvature  parameter 
will  only  be  examined  using  values  found  in  literature. 

Two  figures  of  merit  are  used  to  assess  the  quality  of  the  correlation  in  representing  the  data.  First,  the  average 
percent  difference  (a.p.d.)  is  defined  as  the  difference  between  the  experimental  and  predicted  values  as  a  percentage 
of  the  predicted  correlation  value.  These  are  averaged  over  the  number  of  data  points.  Secondly,  the  number  of  data 
points  (on  a  percent  basis)  that  fall  within  +10%  of  the  predicted  values  indicates  the  degree  to  which  the  correlation 
captures  the  data  trend. 

The  analysis  begins  with  the  general  form  of  the  Nusselt  correlations.  It  applies  the  least-squares  best-fit  method 
to  curved  experimental  data  at  the  inlet  location  (0  =  9°)  for  AR  =  4  across  all  velocities,  heat  fluxes,  and  inlet 
temperatures  while  employing  the  hydraulic  diameter  and  bulk  temperature.  The  resulting  correlation, 

Nu  =  0.118Re  68  Pr26  (4) 

captures  96%  of  the  data  within  +10%  of  the  predicted  values  with  an  average  percent  difference  (a.p.d.)  of  3.3. 
This  correlation  was  generated  over  84  steady  state  data  points.  Similar  results  are  achieved  when  the  analysis  is 
performed  for  the  mid-channel  (6  =  90°)  and  exit  (6  =  171°)  locations,  with  98  %  and  100  %  of  data  points  falling 
within  +10%  of  the  respective  correlations.  This  indicates  that  at  any  one  angular  location  the  data  trend  is  very 
consistent  and  describable  with  each  expression  generated. 

With  the  addition  of  the  viscosity  ratio  term  to  the  general  form  of  the  equation,  no  increase  occurs  in  the 
percentage  of  data  captured  for  the  inlet  and  mid-channel  locations  when  using  the  same  set  of  experimental  data  as 
before.  However,  a  2%  decrease  is  found  for  the  exit  location.  Regardless  of  the  decrease,  all  the  correlations 
generated  thus  far  are  of  good  quality. 
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A  similar  process  is  applied  to  AR  =10  curved  channel  experimental  data  under  the  same  conditions.  Resulting 
correlations  from  the  general  form  at  the  inlet,  mid-channel,  and  exit  locations  captures  95%  or  greater  of  the  data 
within  ±10%  of  the  predicted  values,  respectively.  Corresponding  a.p.d.  values  are  3.6  or  lower.  Figure  3  is  an 
example  of  the  correlation  generated  in  the  general  form  using  mid-channel  data.  These  correlations  are  each 
generated  over  78  steady  state  data  points.  Expressions  with  the  appended  viscosity  term  result  in  approximately  the 
same  percentages  and  average  percent  differences. 


Figure  3.  Nusselt  Number  Correlation  of  Present  Data 

The  next  step  is  to  create  correlations  applicable  for  a  broader  range  of  data  and  assess  their  quality.  This  data 
set  encompasses  all  angular  locations  (6  =  9°,  90°,  and  171°),  heat  fluxes,  and  inlet  temperatures  for  AR  =  4  curved 
channel  tests.  Hydraulic  diameter  and  bulk  temperature  are  used  and  the  data  set  contains  252  data  points.  Both 
expressions  shown  here, 


Nu=0.071Re71Pr'31 


Nu=0.071Re71Pr31 


.0172 


(5) 

(6) 


where  the  percent  captured  decreases  by  50%  and  the  a.p.d.  doubles  for  both  expressions.  The  decrease  in  the 
percent  captured  is  attributed  to  the  initial  development  of  the  thermal  boundary  layer.  The  correlations  created 
from  an  analogous  set  of  234  data  points  for  AR  =10  only  decrease  about  10%.  For  aspect  ratio  10,  the  thermal 
boundary  layer  will  be  smaller  thus  a  smaller  reduction  in  the  percent  captured  compared  to  aspect  ratio  4.  Since  the 
percentages  for  both  aspect  ratios  decrease  this  implies  that  a  term  is  needed  to  address  thermal  development. 

Adding  an  axial  location  term  to  Eq.  6  and  reconsidering  all  angular  locations  for  AR  =  4  to  determine  the  best- fit 
correlation  yields 


(  \ 

.03 

Nu  =  0.066Re71  Pr3 

fh 

10 


V. 0703 


x/D, 


hj 


(7) 


where  the  average  percent  difference  is  10.5  and  52%  of  the  data  is  captured  within  ±10%.  The  addition  of  the  axial 
location  term  does  not  improve  correlation  quality  for  aspect  ratio  4  over  the  results  from  Equations  5  and  6.  This 
indicates  another  term  is  needed  in  the  correlation  to  improve  quality  or  use  only  data  from  the  downstream 
locations  that  is  more  thermally  developed.  This  paper  will  address  the  later  of  the  two  options  in  the  next  paragraph 
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and  the  first  option  later  with  the  addition  of  published  curvature  terms.  As  for  AR  =10,  the  resulting  correlation 
with  its  corresponding  data  set  provides  an  a.p.d.  of  5.2  and  captures  92%  of  the  data  within  ±10%.  The  high 
percentage  for  aspect  ratio  10  shows  that  the  development  of  the  thermal  boundary  layer  has  less  of  an  impact  on 
correlation  quality  when  compared  to  aspect  ratio  4. 

In  order  to  further  increase  correlation  quality,  data  from  the  inlet  location  are  omitted  to  reduce  entrance  effects. 
This  results  in  a  percent  captured  of  99%  with  an  a.p.d.  equal  to  1.9  for  AR  =  4  and  96%  and  3.4  for  AR  =  10.  The 
corresponding  correlations  for  aspect  ratios  4  (168  data  points)  and  10  (148  data  points)  are  shown  below  in 
Equations  8  and  9.  These  are  more  accurate  predictors  of  Nusselt  number,  although  with  the  added  restriction  that 
x/Dh  >15.  Also,  notice  the  exponents  for  the  viscosity  ratio  and  axial  location  terms  have  different  signs.  The 
exponents  for  these  terms  tend  to  be  negative  for  all  expressions  for  aspect  ratio  10  regardless  of  the  data  sets  used. 


Nu  =0.014Re{ 


Pr 


36 


r  Pb  ^ 


1  + 


10 


v 1.42 


x/D, 


hj 


AR  -  4 
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(9) 


Heat  transfer  data  collected  from  straight  channel  testing  have  been  put  through  the  same  correlation  analysis  as 
described  above  for  the  curved  channel  heat  transfer  data.  Correlations  follow  similar  trends  in  the  curved  channel 
correlations,  except  the  addition  of  the  viscosity  ratio  term  to  the  general  form  (Eq.  2)  increases  the  percent  captured 
considerably.  Detailed  straight  channel  analyses  can  be  found  in  Nathman  and  Sturgis  (2008)  and  Sturgis  (2008), 
which  include  other  comparisons  such  as  varying  the  length  scale  and  the  temperature  for  which  to  evaluate  the  fluid 
properties.  Their  results  show  that  employing  the  width  as  the  length  scale  resulted  in  slightly  better  correlations 
with  the  hydraulic  diameter,  but  not  significantly  enough  to  dissuade  use  of  the  well-accepted  hydraulic  diameter. 
Also,  no  difference  in  the  correlation  quality  occurred  when  evaluating  fluid  properties  with  the  film  temperature  as 
opposed  to  bulk  temperature.  These  comparisons  are  not  included  in  this  preliminary  curved  channel  analysis.  This 
analysis  continues  to  apply  the  hydraulic  diameter  and  the  bulk  temperature;  however,  future  analyses  may  revisit 
these  comparisons. 

The  second  aspect  of  the  correlation  analysis  involves  examining  the  effect  of  curvature  terms  on  correlation 
quality.  The  curvature  terms  of  interest  found  in  literature  are  summarized  in  Table  1  and  will  be  inserted  into 
Equation  10.  Also  Eq.  10  will  be  applied  to  each  angular  location  and  downstream  locations  (i.e.  mid-channel  and 
exit  locations).  Sturgis  and  Mudawar  (1999)  and  Sturgis  and  Nathman  (2008)  provide  lists  of  curvature  terms. 
When  affixing  Ito’s  curvature  term  to  the  expression  and  setting  the  value  of  exponent  u  to  1,  such  that 


Nu  =  CRem  Prn 


f  V 

Mb 


<Pc 


(10) 


the  resulting  correlations  capture  above  94%  of  the  data  for  each  angular  location  regardless  of  aspect  ratio. 
However,  applying  Eq.  10  to  just  the  downstream  locations,  the  percentage  drops  to  14%  for  aspect  ratio  4.  For  AR 
=  10  the  percentage  stayed  the  same  (~  94%).  Modifications  to  Eq.  10  such  as  adding  the  thermal  flow  development 
parameter  when  using  larger  data  sets  do  not  yield  better  results. 

Several  authors  found  Ito’s  curvature  term  suitable  for  their  heat  transfer  correlations  for  curved  flow.  Seban  and 
McLaughlin  (1963)  use  the  Dittus  and  Boelter  equation  with  the  film  temperature  and  Ito’s  enhancement  term. 
Performance  is  low  evaluating  Eq.  10  using  the  appropriate  properties.  Also,  applying  their  entire  expression 
captures  below  5%  of  the  present  data.  Rogers  and  Mayhew  (1964)  used  the  same  equation  as  Seban  and 
McLaughlin,  but  instead  used  the  bulk  temperature  which  captured  less  than  15%  of  the  present  data.  Though  both 
of  these  studies  used  water,  the  configuration  was  such  that  the  tube  was  heated  around  the  entire  periphery. 
Taylor’s  (1968)  equation  was  only  applied  to  the  current  data  sets  for  each  aspect  ratio  using  downstream  locations 
because  it  includes  an  axial  location  term.  The  percentage  of  data  captured  was  below  6%  because  the  slope  of  the 
data  and  the  correlation  are  different.  Taylor’s  equation  tends  to  over  predict  the  Nusselt  values. 
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Table  1.  Published  Curvature  Terms 

Author 

Curvature  term,  (pcur 

Comments 

Ito  (1959) 

m  =  Re05 

t  cur 

V1 

D  | 

Ac) 

Limited  to  Re(D/dc)2>6 

Sturgis  and 
Mudawar  (1999) 

0  .046 

(P  —  Re  i  \ 

t  cur  D 

s  \ 

Dh 

2Ro, 

Developed  for  the  exit 
location  and  AR  =  2 

Niino  et  al.  (1982) 

/  \  .04 

f  n  \ 

ro  -Re  -02  h 

( 

1  +  sin  n 

Y  l  \ 

V  | 

Developed  for  rectangular, 
asymmetrically  heated 
passages 

t  cur  ~  l)  , 

\dc) 

Lc+l5Dh)_ 

Pratt  (1947) 

(D  =l  +  3.z 

r  cur 

{-] 

Ac  y 

Developed  for  full-periphery 
heated  coiled  tubes 

Sturgis  and  Mudawar’ s  (1999)  curvature  term  is  inserted  into  Equation  10.  High  percentages  of  data  (~  94%)  fall 
within  10%  of  the  resulting  expression  when  applying  it  to  each  angular  location  and  aspect  ratio.  Percentages  drop 
to  40%  for  AR  =  4  when  examining  data  for  downstream  locations.  Sturgis  and  Mudawar  developed  this  curvature 
term  to  apply  to  the  exit  location,  so  the  inclusion  of  data  at  the  mid-channel  and  exit  locations  is  outside  the 
applicability  of  the  term.  The  percentages  for  AR  =10  remain  high;  approximately  94%  for  the  mid-channel  and 
exit  locations  combined. 


Secondly,  the  correlations  developed  by  Sturgis  and  Mudawar  for  each  angular  location  are  applied  to  the  current 
data  set  at  each  corresponding  location.  None  of  the  correlations  captured  any  data  because  they  over  predict  the 
current  data.  This  was  an  unusual  finding  because  Sturgis  and  Mudawar  conducted  similar  testing  for  aspect  ratio  2, 
but  with  a  different  working  fluid.  For  example,  the  correlation  Sturgis  and  Mudawar  developed  for  the  exit 
location  is  given  in  Eq.  11.  As  an  exercise,  the  values  for  the  coefficient  and  Reynolds  number  and  Prandtl  number 
exponents  are  calculated  from  the  best-fit  method  while  the  curvature  term  exponent  is  fixed  at  unity  for  AR  =  4. 
Equation  12  is  the  resulting  expression  using  the  exit  channel  data  with  bulk  temperature  and  hydraulic  diameter. 
Notice  the  coefficient  and  Reynolds  number  exponent  decrease  causing  the  percent  captured  to  increase  to  100% 
between  Eq.  1 1  and  12.  (These  equations  differ  from  Eq.  10  because  the  viscosity  ratio  term  is  omitted.) 


Nu  =  0. 0302Re808 Pr4 
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Next,  the  curvature  term  from  Niino  et  al.  (1982)  is  evaluated  in  a  similar  manner.  Inserting  their  enhancement 
term  into  Eq.  10  captures  over  94%  of  the  data  when  considering  each  angular  location.  The  results  differ  from 
Sturgis  and  Mudawar  where  98%  of  the  data  is  captured  for  the  mid-channel  and  exit  locations  for  aspect  ratio  4  as 
opposed  to  aspect  ratio  10  where  only  5%  of  the  data  is  captured.  Lastly,  the  full  equation  developed  by  Niino  et  al. 
does  not  capture  any  of  the  current  data  regardless  of  how  the  data  is  sorted. 

Results  similar  to  those  previously  described  are  obtained  when  employing  Pratt’s  (1947)  curvature  enhancement 
term.  Equation  10  with  Pratt’s  curvature  term  produces  the  same  results  as  Ito’s  term  applied  to  data  sets  of  each 
angular  location  and  downstream  locations  from  different  locations  for  each  aspect  ratio. 

Another  way  to  evaluate  the  curvature  terms  from  Table  1  is  graphically.  These  curvature  terms  are  plotted  in 
Figure  4  along  with  curved-to-straight  channel  Nusselt  number  data  for  AR  =  4  at  the  exit  location.  Notice  the 
slopes  are  approximately  the  same  between  the  data  and  the  curvature  terms.  However,  the  published  curvature 
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terms  over  predict  the  data.  Pratt’s  term  does  not  depend  on  Reynolds  number  producing  in  a  horizontal  line.  More 
analysis  should  be  done  by  adding  more  data  from  other  channel  locations  in  order  to  examine  trends  over  a  wider 
range  of  data. 
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Figure  4.  Published  Curvature  Terms 


B.  Pressure  Gradient 

Curvature  effects  on  pressure  gradient  are  examined  experimentally.  Tests  with  no  heat  input  were  performed  for 
AR  =  4  and  10  in  both  curved  and  straight  channels  where  inlet  pressure  was  held  constant  at  2.1  bars  (30  psia). 
Two  pressure  taps  located  near  the  inlet  and  two  located  near  the  exit  (see  Figs.  1  &  2)  measured  the  local  radial 
pressure  gradients.  The  straight  channel  has  four  taps  in  corresponding  locations  to  the  curved  channel. 

Forcing  a  flow  to  bend  or  turn  induces  secondary  flows  which  mean  fluid  velocity  towards  the  concave  wall 
increases  and  fluid  velocity  near  the  convex  wall  decreases  (Kakac  1987).  For  HARCC,  the  concave  wall 
corresponds  to  the  bottom  of  the  flow  area  and  the  convex  wall  is  the  top  of  the  channel.  Curvature  results  in  radial 
pressure  gradients  where  high  pressure  is  located  near  the  bottom  of  the  channel  and  low  pressure  near  the  top.  The 
near-inlet  and  near-outlet  radial  pressure  gradients  are  plotted  in  Figure  5  as  a  function  of  mean  velocity.  Notice  the 
radial  pressure  gradient  strength  increases  as  flow  moves  through  the  curved  channel  for  aspect  ratio  10.  This 
indicates  that  the  secondary  flows  are  becoming  more  developed  and  stronger  which  would  produce  a  larger 
pressure  gradient  at  the  exit. 


The  strength  of  the  pressure  gradient  also  relates  to  centripetal  acceleration.  Figure  5  shows  radial  pressure 
gradients  increase  with  mean  velocity.  This  trend  can  also  be  found  from  a  reduction  of  the  Navier-Stokes  equations 
for  fully  developed  flow  in  a  circular  curved  channel.  With  no  mean  radial  component  of  velocity,  the  Navier- 
Stokes  equation  for  the  radial  direction  in  cylindrical  coordinates  becomes, 


pv2e  _  dP 
r  dr 


(13) 


where  p  is  the  fluid  density,  ve  is  the  tangential  velocity.  A  finite  radial  pressure  gradient  exists  to  balance  the 
centripetal  acceleration  of  the  fluid  which  behaves  like  U2/R0.  The  positive  dP/dr  indicates  that  pressure  increases 
away  from  the  center  of  curvature.  Equation  13  is  shown  in  Fig.  5  using  the  current  conditions  for  AR  =10  and 
substitutes  ve  =  U  and  r  =  Ro.  It  predicts  the  pressure  gradients  near  the  exit  of  the  curved  channel  over  the  velocity 
range  tested. 
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Figure  5.  Radial  Pressure  Gradients  with  Mean  Velocity 

Another  aspect  of  curvature  effects  entails  examining  the  pressure  gradients  in  the  flow  just  upstream  of  the  inlet 
and  immediately  downstream  of  the  exit.  Since  the  flow  is  incompressible,  it  accounts  for  flow  curvature 
downstream  by  forming  a  radial  pressure  gradient  before  the  start  of  curvature.  This  can  be  seen  in  Figure  6,  which 
is  a  CFD-generated  pressure  distribution  for  aspect  ratio  10  at  a  velocity  of  15  m/s  with  an  adiabatic  boundary 
condition  and  is  the  equivalent  to  an  unheated  test.  Flow  is  from  left  to  right  in  the  lower  branch  and  the  curved 
section  begins  at  v  =  0.  The  CFD  results  determine  that  the  radial  pressure  gradient  begins  forming  approximately 
1 1  hydraulic  diameters  upstream  of  the  inlet  and  persists  until  approximately  1 8  hydraulic  diameters  downstream  of 
the  curved  channel  exit  for  the  AR  =  10,  U  =  15  m/s,  and  Tin  =  20  °C  case. 
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Figure  6.  Pressure  distribution  for  U  =  15m/s,  Tin  =  20  °C  with  an 
adiabatic  boundary  condition. 
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IV.  Conclusion 

High  aspect  ratio  curved  and  straight  channel  testing  is  conducted  in  a  test  rig  which  models  a  channel  in  a 
milled-liner  chamber.  Nusselt  number  correlations  are  generated  from  current  AR  =  4  and  10  curved  channel  data  as 
functions  of  Reynolds  number,  Prandtl  number,  viscosity  ratio,  axial  location,  and  curvature.  Correlations  at  each 
angular  location  capture  the  data  trend  well.  Aspect  ratio  4  heat  transfer  correlation  that  includes  all  angular 
locations  predicts  only  52%  of  the  data  within  10%,  where  as  correlation  for  aspect  ratio  10  captures  92%.  Low 
percentages  indicate  thermal  boundary  layer  development.  Improvement  in  correlation  quality  is  found  by  omitting 
inlet  data.  Equations  8  and  9  are  recommended  for  aspect  ratios  4  and  10  respectively  with  the  added  limitation  that 
x/Dh  >  15. 

Several  published  curvature  terms  are  appended  to  the  correlations  generated  with  curved  channel  data.  With  the 
addition  of  most  curvature  terms,  the  correlations  predicted  over  90%  of  the  data  at  each  angular  location.  However, 
when  an  entire  published  correlation  was  applied  to  the  data,  the  data  trends  were  not  well  captured. 

Curvature  effects  on  pressure  gradients  are  examined.  The  strength  of  radial  pressure  gradients  directly  relates  to 
the  centripetal  acceleration.  A  reduced  form  of  the  Navier-Stokes  equations  predicted  the  near  exit  radial  pressure 
gradient  well.  CFD  pressure  distributions  showed  the  effects  of  curvature  upstream  and  downstream  of  the  curved 
channel.  Radial  pressure  gradients  begin  forming  approximately  1 1  hydraulic  diameters  upstream  of  the  inlet  and 
persist  until  approximately  18  hydraulic  diameters  downstream  of  the  curved  channel  exit. 

Future  work  entails  completing  this  series  of  testing  (AR  =1,2,  and  8)  and  modifying  the  channel  to  incorporate 
more  instrumentation.  CFD  work  is  currently  underway  to  investigate  the  formation  of  secondary  flows  in  HARCC. 
Also,  a  program  is  in  place  that  will  test  at  rocket  like  conditions. 
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Purpose 


•  Develop  design  tools  for  use  with  high  aspect  ratio 
coolant  channels  (HARCC) 

•  Contribute  to  the  thermal  management  knowledge  base 
of  regeneratively  cooled  liquid  rocket  engines  by 

-  Developing  correlations 

-  Identifying  trends 

-  Analyzing  the  fin  effects  of  the  channel  walls 

-  CFD  for  visualization,  understanding,  and  insight  into  flow 

•  Create  database  with  straight  and  curved  channel  data 
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Strategy 


(Current  effort) 


•  Phase  I  -  Water  Flow  Loop  Testing 

-  Straight  and  curved  channels 

-  Rich  experimental  database  (AR,  Re,  q”,  Tin,  curvature) 

-  Develop  correlations,  predictive  tools 

-  CFD  for  insight  into  flow 

-  Limitations:  low  Re  and  q”  for  rocket  applications 


Current 
^  effort 


j 


•  Phase  II  -  Hot-Fire  Testing 

-  Sub-scale  combustion  chamber  with  groups  of  different  sized  coolant 
channels 

-  Achieve  rocket-like  q”  and  Re,  multiple  coolants,  validate  tools 
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Coolant  Channel  Model 


T- insert 


•  Test  channel  models  a 
rectangular,  one-side-heated 
regen  channel  with  conducting 
side  walls 


channel 
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coolant 
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tt  t  m  1 1  tt 
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Milled-liner  combustion 
chamber  cross-section 


insulation 


copper 
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heater 


insulation 


Test  channel  cross-section 
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Test  Channel  Design 


T/Cs  in  neck  measure  temp  gradient,  hence  q” 
T/Cs  embedded  in  side  walls  (fins) 

Heat  input  adjusted  by  voltage  to  heaters 
Channel  width  to  land  width,  W:tf  =  1 :1 
One-sided  heating,  conductive  side  walls 


T-insert 
(k  =  5  W/m-K) 

Copper  channel 
(k  =  390  W/m-K) 


y2  tf 


T/Cs' 


W 


AR=4 

AR=2 

AR=1 


cartridge 

heater 
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Instrumentation 


Length  =  ttR0  =  15  in 
U  =  5,  10,  15  m/s 
Re  <  300,000  {Rew<  180,000) 
Tin  =  20  or  35  C 
4.3  <  Pr  <  7 
AR  =  4  and  10 
Curved  Channel 

q”  <  0.8  Btu/in2-s  (130  W/cm2) 
21  <  2  Ro/Dh<  24 
6,000  <  De<  66,000 
21  <  g*<  190 


Copper 

channel 


Flow  area 


Pressure 

taps 


External 
support  rib 


Straight  Channel 
•  q”  <1.1  Btu/in2-s  (180  W/cm2) 


Therm  ocouples 


DISTRI 


Cart  ridge 
heaters 


Curved  Test  Channel  in  Flow  Loop 

Pressure  Cartridge 
Taps  Heater  Slots 


Aluminum  Copper  Thermocouples 

Entry  Length  Test  Channel 
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Straight  Test  Channel  in  Flow  Loop 


Pressure  Transducers 


Inlet  Aluminum  Copper  Cartridge  Exit 

Plenum  Entry  Length  Test  Channel  Heaters  Plenum 


DISTRIBUTION  A,  Approved  for  public  release. 


9 


Data  Reduction 


•  Local  average  convection  coefficient,  h , 
determined  over  3  walls  at  each  axial  location 
along  channel 

•  Convection  coefficient  calculated  using  an  energy 
balance 

Heat  input  =  Heat  convected  to  coolant  +  Heat  loss 

•  Used  in  calculating  Nusselt  numbers  (Nu): 


k 
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Correlation  Analysis 


•  Goal:  Generate  a  design  and  analysis  tool  that 
accurately  predicts  the  heat  transfer  coefficient  in  a 
coolant  channel  (i.e.  a  Nusselt  Number  Correlation) 


Nu=  C  Rem  Prn 
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Mb 

\Pw  J 
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1  + 
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cur 
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•  Correlate  Nu  using  a  least-squares  best  fit 
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Correlation  Deduction 


•  Data  sets  include: 

-  Each  aspect  ratio  (4  or  1 0) 

-  All  heat  fluxes,  velocities,  and  inlet  temperatures 

-  3  axial  locations  (inlet,  mid-channel,  exit,  and  combinations  of 
the  3) 


•  Parametric  Analysis  of  Nusselt  correlations: 

-  Addition  of  terms  (example:  viscosity  ratio) 


•  Assess  the  quality  of  the  correlations 

-  Average  percent  difference  (a.p.d.) 

-  Percentage  of  data  points  that  fall  within  +/- 10%  of  the  predicted 
values  (%  captured) 
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Correlation  Deduction  (cont.):  AR  =  4 

•  Conditions: 

-  Curved  channel  data,  hydraulic  diameter,  bulk  temperature,  all  inlet  temps,  AR  =  4 

•  Each  individual  angular  location 

-  Correlations  capture  data  well  (~  96-100%  captured)  (84  data  points) 

•  All  angular  locations  (Inlet,  mid-channel,  and  exit) 

-  Correlations  capture  ~  50  %  (252  data  points) 

-  Thermal  Boundary  Layer  Development  -  omit  inlet  data 

•  Only  downstream  locations  (mid-channel  and  exit  locations,  x/Dh  >15)  (168  data  points) 
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orrelation  Deduction  (cont.):  AR  =  1 

onditions: 

-  Curved  channel  data  , hydraulic  diameter,  bulk  temperature,  all  inlet  temps, 

•  Each  individual  axial  location 

-  Correlations  capture  data  well  (~  96  %  captured)  (78  data  points) 

•  All  axial  locations  (Inlet,  mid-channel,  and  exit) 

-  Correlations  capture  ~  89%  (234  data  points) 

-  Thermal  Boundary  Layer  Development  -  omit  inlet  data 

•  Only  downstream  locations  (mid-channel  and  exit  locations,  x/Dh  >15)  (148  data  points) 

(  a.p.d  /  %  Captured  ) 
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[relation  Deduction  (cont.):  AR=  10 
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Straight  Curved 

(430  data  points)  (1 68  data  points) 


Straight  vs.  Curved 


•  Conditions:  hydraulic  diameter,  bulk  temperature,  all  inlet  temps,  AR  =  4 

•  Only  downstream  locations  (mid-channel  and  exit  locations,  x/Dh  >  15) 
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Parametric  Results 


•  Addition  of  terms 

-  Axial  location  term  increases  correlation  quality 

•  Recommended  Correlations  (mid-channel  and  exit  locations,  x/Dh  >  15) 
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•  Curved  vs.  Straight 

-  Different  slopes 

-  Viscosity  ratio  increases  correlation  quality  more  for  straight  channel  data 

•  Correlation  Limits 
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Addition  of  Curvature  Terms 


•  Appending  published  curvature  terms  did  well  at  each  location 

•  Low  performance  for  all  locations 


Nu  =  C  Rem  Prn 
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Developed  for  full-periphery 
heated  coiled  tubes 
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Curvature  Terms  (cont.) 


•  Addition  of  curvature  terms  at  each  angular  location  capture  data 
trends  well  for  AR  4  and  1 0  (~94%) 

•  Percentages  low  when  evaluated  using  all  angular  locations 

•  Ito’s  term: 

-  1 4%  for  mid-channel  and  exit  angular  locations  for  AR  4 

-  94%  for  mid  and  exit  angular  locations  for  AR  1 0 

•  Using  Pratt  term  had  same  results  as  Ito 
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Curvature  Terms  (cont.) 


•  Sturgis  and  Mudawar 

-  40%  for  downstream  angular  locations  for  AR  4 

-  94%  for  mid  and  exit  angular  locations  for  AR  1 0 

•  Niino 

-  98%  for  mid  and  exit  angular  locations  for  AR  4 

-  5%  for  mid  and  exit  angular  locations  for  AR  1 0 
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Published  Curvature  Terms 
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Curvature  Effects: 
Radial  Pressure  Gradients 
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♦  Near  Inlet 

■  Near  Exit 

■  —  Navier-Stokes  Eq. 
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fully  developed 
flow  in  a  circular 
curved  channel 

No  mean  radial 
component  of 
velocity 

Navier-Stokes 
equation  for  the 
radial  direction 
becomes: 


pv2e  _  dP 
r  dr 


DISTRIBUTION  A,  Approved  for  public  release. 


22 


(uiui)a 


Curvature  Effects: 

Radial  Pressure  Gradients  (cont.) 


■ 

■ 

P(psia):  14  18  22  26  30  34  38  42  46  50 

X(mm) 


Conditions: 
AR=  10 
U=  15  m/s 
Tin  =  20  QC  case 


•  Radial  pressure  gradient  begins  forming  ~  11  Dh  upstream  of  inlet 

•  Persists  ~  1 8  Dh  downstream  of  the  exit 
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Summary 


i#  . wj 

Correlation  Analysis: 

-  Correlations  generated  at  each  location  captured  data  trend  well 

-  Correlations  with  viscosity  ratio,  and  x/D  terms  capture  over  96%  of  data 
within  1 0%  of  the  prediction  (x/Dh  >  1 5) 

-  Addition  of  published  curvature  terms  captured  data  trend  well  at  each 
location 

-  Correlation  able  to  predict  heat  transfer  using  design  inputs,  however 
limited  to  current  test  conditions 


Curvature  Effects: 

-  strength  of  radial  pressure  gradients  directly  relates  to  the  centripetal 
acceleration 

-  A  reduced  form  of  the  Navier-Stokes  equations  predicted  the  near  exit 
radial  pressure  gradient  well 

-  CFD  radial  pressure  distributions  showed  the  effects  of  curvature 
upstream  and  downstream  of  the  curved  channel 

•  Radial  pressure  gradients  begin  forming  ~  1 1  Dh  upstream  of  the  inlet 

•  Persist  ~  18  hydraulic  diameters  downstream  of  the  exit  (AR  10,  15  m/s,  20QC) 
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Future  Plans 


•  Curved  Channel  Testing 

-  Complete  AR  =  1 , 2,  and  8 

-  Revisit  comparing  to  straight  channel  data  -  Curvature 
Enhancement 

-  Assess  friction  factor 


•  Continue  CFD  effort 


•  Extend  database  to  higher  Re,  higher  q”,  different  fluids 
(Phase  II) 
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Questions  ? 


Data  Reduction 


Convection  Coefficient 


q" (W + tf)  =  Xh  (Tj  -  Tb)  As,  +  q'flntCond 

\  7  \  \  \ 


3  T/Cs 

in  neck  average  h 

for  cross-section 
(3  walls) 


bulk  temp 
assumption 

fin  temperature  data 


estimated  from 
fin  temp  gradient 


# 
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